Multiple Ionization Bursts in Laser-Driven Hydrogen Molecular Ion 
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Theoretical study on H J in an intense infrared laser field on the attosecond time-scale reveals that 
the molecular ion shows multiple bursts of ionization within a half-cycle of the laser field oscillation, 
in contrast to the widely accepted tunnel ionization picture for an atom. These bursts are found to 
be induced by transient localization of the electron at one of the nuclei, and a relation between the 
time instants of the localization and the vector potential of the laser light is derived. Furthermore, 
an experimental scheme is proposed to probe the localization dynamics by an extreme ultraviolet 
laser pulse. 

PACS numbers: 32.80.Rm, 33.80.Rv 



Many interesting and useful phenomena induced by the 
intense laser-matter interaction, such as higher-order har- 
monics and attosecond pulse generation as well as non- 
sequential double ionization, are known to be initiated 
by the release of an electron from the parent atom or 
molecule [l[. This ionization process in an intense laser 
field is often understood in terms of the quasi-static tun- 
nel ionization picture. According to this picture, the elec- 
tron tunnels through the barrier created by the combi- 
nation of the binding potential of the ionic core and the 
electric potential of the laser field. Thus, in the oscillat- 
ing electric field of laser light, the electron is expected 
to escape most likely at the peak of the electric field 
strength in every half-cycle when the barrier becomes 
the thinnest. This expectation is in agreement with the 
temporal dependence of the ionization rate for an atom 
predicted by the well-known tunneling formula 0, H[ and 
shown in Fig. [TJa). In a recent experiment, the ioniza- 
tion yield of Ne + was measured as a function of time 
with attosecond resolution, and the results also support 
the picture above 0]. However, results of our present the- 
oretical study indicate that this popular picture needs to 
be modified for molecules: they can exhibit not only a 
single but multiple bursts of ionization within a half-cycle 
of the laser field. For example, the ionization rate of 
shown in Fig. QJb) has two maxima within a half-cycle 
and minima near the peaks of the laser electric field. 

In this Letter we show that these multiple ionization 
bursts (MIBs) are related to transient electron localiza- 
tion at one of the protons on the attosecond time scale. 
The localization is due to the strong coupling between a 
pair of states with opposite parities, called the charge- 
resonant (CR) states Q, which leads to more complex 
electron dynamics in a molecule than in an atom, as it 
has been reported before @, Q • Despite the complexity of 
the electron dynamics, we derive a simple relation of the 
instants of electron localization to the vector potential of 
the laser field and the coupling strength between the CR 
states. This enables us to generalize our observations for 
the simplest molecule and to predict a similar modifica- 
tion of the tunnel ionization picture for other molecules 
having the CR states as well. To probe the transient elec- 
tron localization in experiments, we propose a method to 
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FIG. 1: Time-dependent ionization rates (solid blue curves) 
of (a) H atom 0] and (b) molecular ion over a single period 
at the peak of a laser pulse, represented by its electric field 
(dashed red curves). The ionization rate of is for the range 
of internuclear distance, 6.75 a.u. < R < 7.25 a.u., and was 
obtained by numerical simulation of TDSE. Pulse parameters 
in the calculations were wavelength: 800 nm, peak intensity: 
6 x 10 13 W/cm 2 , pulse duration: 26.69 fs (FWHM). 



retrieve the laser-dressed quantum state of from an 
interference pattern in the photoelectron momentum dis- 
tribution generated by applying an additional attosecond 
extreme ultraviolet (XUV) pulse. 

For the present analysis wc have solved the time- 
dependent Schrodinger equation (TDSE) for a model of 
in which the electronic and nuclear motions are re- 
stricted along the polarization direction of the linearly 
polarized laser pulse. The Hamiltonian of this system is 
given by (Hartree atomic units are used throughout) [U : 



H(t) = 



_i d?_ j_ 

dR? ~ 2\i z dz 2 
1 1 



(1) 



■PzE(t), 



where R is the internuclear distance, z is the position 
of the electron measured from the center-of-mass of the 
protons, fi R = M/2 and [i z = 2M/(2M + 1) are the re- 
duced masses with M = 1.836 x 10 3 a.u. being the proton 
mass, a e = 1.0 a.u. and a n = 0.03 a.u. are the soft-core 
parameters 0, /3 = 1 + 1/(2M+ 1), and E(t) is the laser 
electric field. It has been shown before that this type of 
reduced-dimensional models reproduce experimental re- 
sults at least qualitatively jfj|-KLl| . The ionization rate in 
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Fig.QJb) was calculated as the out-going probability flux 
from the rectangular domain, —8.4 a.u. < z < 8.4 a.u. 
and 6.75 a.u. < R < 7.25 a.u., toward \z\ oo nor- 
malized by the probability inside this domain, under a 
laser pulse with the wavelength of 800 nm, peak inten- 
sity 6 x 10 13 W/cm 2 , and FWHM duration 26.69 fs. 

We present in Fig. [5Ja) how the electron density 
evolves in time under the same laser pulse and in the 
same scope of R as used in Fig. [IJb) . This result shows 
that the electron density bound to the protons at around 
z « ±3.5 a.u. is released in several bunches within a 
half-cycle of the laser field, consistently with the obser- 
vation in Fig.[TJb). In the following, we trace the origin of 
these MIBs in the time-evolution of the electron density 
by simplifying the model. 

Figure HJh) shows the electron density calculated by 
fixing the nuclear positions at R = 7 a.u. and letting 
AI — y oo. This fixed- nuclei model reproduces the re- 
sult of the moving-nuclei model in Fig. [2ta) almost per- 
fectly, indicating that the coupling between the elec- 
tronic and nuclear motions is not essential for the for- 
mation of the MIBs. In Fig. [2fc), in the fixed-nuclei 
model, we absorbed the ionizing wavepackets soon after 
they left the protons by using a cos 1//2 -mask function set 
over 6.6 a.u. < \z\ < 11 a.u. This result shows that the 
bound electron density is transiently localized at one of 
the two protons at the same instants as the ionization 
bursts marked by the circles in panel (b). Part of the 
wavepackets released in the bursts during the time pe- 
riod in which \E(t)\ is decreasing (—0.25 < t < and 
0.25 < t < 0.5 laser cycles in Fig. [2]) can return to the 
protons P, [l3] • Such rescattering wavepackets create the 
modulation of the density including the enhancements at 
the cross marks in Fig. [2jb) due to the interference with 
the bound wavepacket and also the finer fringes in the 
ionizing density at \z\ 8 a.u. due to the interference 
with the wavepackets just released from the protons. 

From the analysis so far, we can conclude that the 
MIBs are induced by the transient electron localization 
at one of the protons while the (rescattering) dynamics 
of the electronic wavepackets in the continuum does not 
affect the time instants of the MIBs. This also substan- 
tiates that the present ID models are sufficient for the 
analysis since effects of spreading in the transverse direc- 
tion are negligible for the bound wavepackets. In previ- 
ous theoretical results, the laser-driven electron dynamics 
inside the molecule has been recognized to be complex 
and sometimes counter- intuitive [fj, and, in particular, 
the existence of the sub-laser-cycle electron localization 
has been pointed out 0j. The present work elucidates 
that this ultrafast electron localization also makes the 
ionization dynamics of qualitatively different from 
that of atoms and may give rise to a paradigm change 
from the widely accepted tunneling ionization picture. 

In order to further analyze the phenomenon of MIBs, 
we investigate the bound state dynamics using a simple 
model which incorporates only the ground (|g)) and first- 
excited (|u)) electronic states of H^". These two states 
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2: Electron density (a) integrated over 6.75 a.u. < R < 
a.u. in the moving-nuclei model and (b) in the fixed- 
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nuclei model at R — 7 a.u. as a function of time and electron 
position (shown on a log 10 -scale). (c) Electron density in the 
fixed-nuclei model with absorber near the nuclei, (d) Electron 
populations in the localized states, |L) and |R), in the two- 
state model, (e) Vector potential (solid line) and electric field 
(dashed line) of the laser light. Horizontal and vertical grid 
lines represent A(ti oc ) and t\ oc , respectively. 



can be written as superpositions, |g) = [|L) + |R)]/v2 
and |u) = [|L) — |R)]/\/2, of the ground states of H 
atom, |L) and |R), centered at the protons at z = —R/2 
and z = R/2, respectively. During the interaction with 
laser light, the cigenstates |g) and |u) mix, and the elec- 
tron is driven between the two protons. This generates a 
large transition dipolc d gn = — (g\z\u), which is propor- 
tional to R at large R, and is known as the CR mecha- 
nism (H, |9|, [l3j . The two-state model is applicable to the 
present case since at R = 7 a.u. the |g) and |u) states 
are almost degenerate and well isolated from the higher- 
lying excited electronic states on the scale of the IR pho- 
ton energy. The time-evolution of the state \^(t)) of this 
model in a laser field can be described in the basis of 
Floquet states I14 - fl6l |. Using a series expansion for the 
Floquet states [T(j, we obtained expressions for the lo- 
cal populations at the respective protons, |(L|\l/(f))| 2 and 
|(R|*(t))| 2 . From the condition d[\(L\^(t))\ 2 ]/dt = 0, we 
derive that population at one of the protons is maximized 
at the instants t\ oc at which the vector potential satisfies 

A(t\ oc ) = A sin(wti oc + (fi) = —— , (2) 

where m = 0, =bl, ±2, . . . , and ^ is a mixing an- 
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gle determined by the quantum amplitudes c\ and 
c 2 of the Floquet states which reduces to |g) and 
|u) at zero laser intensity, respectively, with cosx = 
(| ci| 2 - M 2 )/C, sin X = 2 Im[ Cl c 2 ]/C, and C = 
\/(l c i| 2 — l c 2| 2 ) 2 + 4(Im[c*C2]) 2 . Furthermore, from the 
sign of cP\(L\i$f(ti oc )}\/dt 2 , we can predict that the elec- 
tron density is localized at the proton down the slope of 
the electric potential of the laser field at those t\ oc corre- 
sponding to odd m. Previously, the electron localization 
was analyzed in terms of the phase-adiabatic states, but 
numerical computation was necessary to predict the in- 
stants of the maximum localization [7j. 

In Fig. dfd), the local populations, |(L|*(t))| 2 and 
| (R|5'(i)) | 2 , calculated by the two-state model are plotted 
as a function of time. Comparison with Fig. [2jc) shows 
that the two-state model does reproduce the bound elec- 
tron dynamics obtained by the numerical simulation for 
the ID fixed-nuclei model. Therefore, the number and 
time instants of the sub-cycle ionization bursts can be 
traced back to the simple condition ©. 

We may emphasize that the strong and exclusive cou- 
pling between a pair of states, and hence the trapping 
of the population within such a pair, are the essence 
of the counter-intuitive two-state dynamics inducing the 
MIBs. Thus, the current analysis applies in general to 
molecules in which the ground state is strongly coupled 
to one particular excited state, and to laser wavelengths 
which do not coincide with a resonant transition to any 
other state. In atoms, the first-excited state is coupled 
not only to the ground state but also with higher-lying 
excited states. Therefore, the electron localization can 
not be induced, and MIBs have not been observed. Even 
in H^~, at around the equilibrium internuclear distance 
(R cq ~ 2 a.u.), the energy levels of |g) and |u) are not 
well isolated from the other states, and the transition 
dipole is not exclusive between these two states [!, IH, [l3[ . 
Therefore, the ionization dynamics becomes similar to 
atoms, and the standard quasi-static tunnel ionization 
picture should be recovered. On the other hand, in 
at larger R, where the coupling strength d gu ~ R/2 be- 
comes large, the electron density can be localized more 
than twice within a half-cycle at rather moderate laser 
intensity. Note that the coupling between CR states can 
be large and exclusive already at the equilibrium struc- 
ture in other molecules [l7| ■ 

As shown above, MIBs are caused by the localization 
of the electron at one of the protons. Previously, theo- 
retical studies have shown that the time-evolving asym- 
metry of the electron density inside can be probed by 
an attosecond XUV pulse via the asymmetry of the pho- 
toelectron yield and momentum in opposite directions 
along the XUV laser polarization parallel to the molecu- 
lar axis [H, [l9| . We propose here to extend this scheme 
by setting the polarization of the XUV pulse at an angle 
to the molecular axis and the polarization of the IR driv- 
ing pulse and by analyzing the 2D photoclcctron momen- 
tum distribution. As we will show below, this potentially 
enables us to reconstruct the temporal evolution of the 



amplitudes ct,(t) and c R (t) of the two quantum states |L) 
and |R) composing the laser-dressed state of . 

To this end, we model the photoionization process due 
to the XUV pulse as a one-photon transition from the IR- 
laser-dressed bound state, \^f(t)) = Cl(£)|L) +cr(£)|R), 
to the Volkov state of drift momentum k in the IR laser 
field. For simplicity, we assume here that the XUV pulse 
is polarized perpendicular to the molecular axis and its 
electric field has a Gaussian envelope as Exuv(i) = 
exuv^xuv exp [-(t - io) 2 / T xuv] cos(ilt + ^xuv)- At 
the limit of Txuv <C 2ir, where rxuv is the pulse du- 
ration of the XUV pulse and Tir is the period of the IR 
driving pulse, by applying Laplace's method of asymp- 
totic analysis for the S- matrix element, we obtained the 
photoclcctron momentum distribution as 



2_2 

XUV 



\S fi (k,t )\ 2 ~ ^[E^t 

^ Hq(MI 2 e, 



x exp 



'XUV 



datom(q(to)) 



(3) 



2 2 

x {| CL (to)| 2 + |cR(to)| 2 

+ 2Re [c L (t )c R (to)] cos[q(t ) ■ R] 
- 21m [cL(*o)c R (*o)]sin[q(t ) -R]}, 

where q(to) = k + A (to) is the photoelectron momentum 
k offset by the vector potential A (to) of the IR field at the 
peak to of the XUV pulse, E g and E u are the energies of 
|g) and |u), respectively, and et atom (q) = (e iq ' r |r|0i s (r)) 
is the transition dipole between H atom Is state to the 
plane wave state. 

By fitting the model eq. ([3]) to an experimentally ob- 
tained photoelectron momentum distribution, the quan- 
tum amplitudes {ch(to), CR,(to)} at the peak of the XUV 
probe pulse can be retrieved. In order to demonstrate 
the accuracy of this retrieval procedure, we simulated 
momentum distributions by solving the TDSE for a 2D 
model of with the internuclear distance fixed at R = 7 
a.u. To fulfill the condition rxuv ^ ?ir, we used an IR 
laser pulse with wavelength of 1400 nm, peak intensity 
1.5 x 10 13 W/cm 2 , and FWHM pulse duration 14.01 fs. 
For the XUV laser field, the peak intensity was set at 
1.0 x 10 12 W/cm 2 , wavelength at 25 nm, and FWHM at 
500.3 as. By changing the delay At from the peak of the 
IR pulse to that of XUV pulse, photoelectron momentum 
distributions were calculated. Then, from these distribu- 
tions, the background momentum distribution calculated 
by applying only the IR laser pulse was subtracted. 

Figure [3] shows the XUV photoionization signal ob- 
tained by the numerical simulations at six delay times 
from At = to 0.25 IR laser cycles. As the Gaussian fac- 
tor in the model formula © suggests, the photoelectron 
momentum is distributed around the energy-conservation 
circle of radius y / 2[(E g + E u )/2 + Q], whose center is 
streaked by the vector potential of the IR laser field as 
k = — A(t ) [20|. This ring distribution is multiplied by 
the atomic p^-wave factor due to |ct a tom| 2 , as well as the 
two-center interference pattern from which the ultrafast 
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FIG. 3: XUV photoelectron signals at the XUV-IR delays of 
(a) At = 0.00 (b) At = 0.05 (c) At = 0.10 (d) At = 0.15 (e) 
At = 0.20 (f) At = 0.25 IR cycles. (R || x and e XU v || y) 
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FIG. 4: (a) Normalized local populations Pl and Pr and (b) 
relative phase o-lr retrieved by the fitting procedure (mark- 
ers) in comparison with the values obtained from the TDSE 
solutions (gray lines). 



evolution of cl(£o) and cr(£o) can be retrieved. 

In actual experiments, it is, however, difficult to de- 
termine the absolute value of the photoelectron momen- 
tum distribution. Under this restriction, the informa- 
tion about the common scale of |cl(£o)| and |cr(£o)| is 
lost, but we can still retrieve the normalized populations, 
Pl = |cl|7(|cl| 2 + |c R | 2 ) and P R = |c R | 2 /(|c L | 2 + |c R | 2 ), 
as well as the relative phase c*lr = arg[cLC R ]. Figure 2] 
shows these quantities retrieved from the distributions 
in FigO (and those at other XUV-IR delays). We also 



computed the exact time-evolution of Pl, -Pr, and qlr in 
the TDSE simulation with only the IR pulse applied, and 
show the results in Fig. @] by gray lines as reference. The 
comparison shows that the present method allows one to 
reconstruct the attosecond time-evolution of the laser- 
dressed quantum state in Hj. We may note that due 
to the present assumption that the polarization direction 
of the XUV pulse is perpendicular to the molecular axis 
we obtain the same momentum distribution when Pl is 
exchanged with Pr. An extension to other alignment an- 
gles is straightforward, and in this case the set of values 
{Pl,Pr,q:lr} can be retrieved uniquely. 

In summary, we have found that there can be MIBs 
from H^" and other molecules with CR states within a 
half-cycle of the laser field oscillation in contrast to the 
widely accepted tunnel ionization picture. These bursts 
have been shown to be induced by the transient electron 
localization inside the molecule on the attosecond time- 
scale, and a simple expression to predict the number and 
instants of the electron localization has been presented. 
The time-evolution of the quantum state exhibiting such 
a localization behavior can be reconstructed in an exper- 
iment with an attosecond XUV pulse. 

Beyond the paradigm change from the tunneling ion- 
ization picture, our findings also suggest that in a HHG 
process, where the intermediate cation is actually dressed 
by the intense driving field, the ultrafast time-evolution 
of the electronic state can be drastically different from 
that under the field-free condition if the cation is in a 
CR state. This should influence the results of techniques 
to reconstruct molecular orbitals [2l[. The MIBs also 
offer a control knob of the attosecond pulse generation 
for which the harmonics at the cut-off of the spectrum is 
used. The photon energy in HHG corresponds to the ki- 
netic energy of the electron acquired in the driving field, 
which is determined by the phase of the field at which the 
electron is released [l|, [l2[ . Our findings suggest that the 
ionization probability at a particular phase, and hence 
the harmonic efficiency at the cut-off, can be controlled. 
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